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ABSTRACT 

We have used a Lagrangian, fully implicit, one-dimensional 
hydrodynamic computer code to evolve thermonuclear runaways in the 
accreted hydrogen rich envelopes of 1.0 M neutron stars with radii 

w 

of 10 km and 20 km. Our simulations produce outbursts which last 
from about 750 seconds to about one week. Peak effective tenoera- 

7 4 

tures and luminosities were 2.5 x 10 K and 8 x 10 L^ for the 10 km 

study and 5.3 x 10“ K end 600 L fcr the 20 km study. Hydrodynamic 

0 

expansion on the 10 km neutron star produced a precursor lasting 
-A 

about 10 seconds. 


INTRODUCTION 

We have studied thermonuclear runaways in the accreted hydro- 
gen rich envelopes of 10 km and 20 km neutron stars using a fully 
implicit, Lagrangian, hydrodynamic computer code which incorporates 

a nuclear reaction network^ We have assumed that the Bursters and 
Transient X-ray sources occur as a result of mass transfer from a 
secondary onto a neutron star in an analogous fashion to the nova 

ohenomena'^ . Reviews of oublished work on this subject can be found 
3 4 

in the literature ’ . 

The published work^ has produced simulations of the Burster 
phenomena which are in reasonable agreement with the observations 
but have not yet reproduced the full range of observed behavior. 

By tr.is we mean that observed Bursts show a very wide range of time 



scales lasting froir. seconds to minutes and the published work has, 
so far, only addressed the characteristics of the shortest bursts. 
If one wants to prolong the time scale, then there are three pos* 
Sibil Itles all of which Involve Increasing the amount of fuel 
available for nuclear burning. First, one can assume that the 
luminosity of the neutron star is low and the mass accretion rate 
Is low and build up a thick hydrogen envelope. Second, one can 
assume a rapid Inflow of material and include the accretion energy 
In the radiative losses; and, third, one can assume that at some 
accretion rate the hydrogen will burn stably and a thick layer of 
helium can be built up on the neutron star. 

3 4 

Because of the success of the studies of Joss and Taam re- 
viewed Ip this volume by Joss, we have concentrated on a regime 
of the (H,L) plane not yet studied by these investigators: that of 
low Internal neutron star luminosity and low mass accretion rates. 
In our case we are trying to both model the long time scale pheno- 
mena observed from some Bursters and to understand the cause of 
the very long time scale outbursts of the Transients. 

Our computer code has previously been used In studies of the 

nova outburst,®’ ’ thermonuclear runaways In the accreted hydro- 
gen rich envelopes of white dwarfs, and is Ideally suited to this 
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Fig. 1. The temperature as a 
function of time. 


Fig. 2. The 
function of 

luminosity as a 
time. 



study. The physics that we include is described elsewhere^. We 
include the p-p chain, the CNO reactions, the triple-a reaction, 

the ^^C(a,Y)^®0 reaction and, finally, assume that the ^^0(a,p)^^F 
15 IQ 

and 0{a,Y) Ne reactions are acting in the fashion described by 

O 

Wallace ano Woosley . We use these last two reactions to provide 
a measure of the rate of depletion of the CNO nuclei during the 
evolutionary sequences and then use this rate of depletion to cal> 
culate an energy generation assuming that a few further proton 
captures will occur. 

RESULTS 

The radius of model 1 was 10 km and it had an envelope mass of 
1.5 X 10“^^ Mg. This is a thick envelope for a neutron star of 
this mass but can be obtained for low accretion rates onto neutron 

f 

stars with low internal temperatures'*. Although some nuclear burn- 
ing could have occurred in the deeper layers during the accretion 
process, we have neglected this and assumed a sharp composition 
interface. We shall refer to this boundary as the core-envelope 
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Fig. 4. The hydrogen (x) and helium Cy) abundances (mass 
fraction) as a function of mass at two times in the evolution 
of the 20 kn neutron star. 
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interface (C£I). We chose a luminosity for the neutron star of 

0.1 which results in an effective temperature of 8 x 10^ K. The 

® 7 S 

temoerature and density at the CEI are 3.3 x 10 K and 3.1 x 10 gm 

•3 

cm respectively. 

We begin the evolution by turning on the nuclear reactions 

and it takes this sequence about 760 seconds of evolution to reach 

g 21 

a temperature of 2.45 x 10 K and an energy generation of 10 erg 

-1 -1 4 

gm s . The surface luminosity has reached 8 x 10 L and the 

7 ^ 

effective tanperature is 2.6 x 10 K(kT '\.2.2 kev). The entire 
envelope has become convective and is mixing fresh unburned nuclei 
into the shell source from the surface. 

Because of our initial conditions, which may be physically 
realizable, the Initial region of peak temperature is not at the 
-12 

CEI but about 10 closer to the surface. The reason for this 

effect is that the electron degenerate conductivity is large 
enough to transport the energy produced in the deepest layers of 
the accreted envelope into the core. This effectively keeps these 
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Fig. 6. The luminosity of the 
neutron star as a function of 
time. 
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regions cooler than the less degenerate regions closer to the sur- 
face. However, the energy produced by the shell source is now 
heating the hydrogen rich layers just below it and their tempera- 
ture slowly climbs until they reach about 3 x 10® K. At this 

9 

point they flash to a temperature of 2.8 x 10 K and an energy 

generation rate of 10 ® erg gm“ s“ . This flash produces an over- 
pressure of a few percent which results In a shock wave that 

reaches the surface some 5 x 10"® sec after the flash occurred. 

The surface luminosity climbs to 2 x 10® L. and the effective 

7 ® 

temperature to 3 x 10 K. The shock causes the surface zones to 

expand at velocities of 2 x 10® km sec~^ . However, this expansion 

-6 

lasts for only about 10 sec and a total excursion in radius of 
1 . 5 meters . 

The shell source slowly moves inward and these layers flash 
to high temperatures although no more shocks occur. The tempera- 
ture as a function of time for the deepest hydrogen rich shell is 
given in Figure 1. The hydrogen burning layers stay hot until all 
of the helium and most of the hydrogen is consumed. It takes 
about 400 seconds for this to occur and during this entire phase 
of nuclear burning the envelope has slowly exoanded to about 12 km 
and remained at virtually constant luminosity. The light curve 
for this simulation is given in Figure 2. 

Once the fuel has been burnt, the radius of the envelope 
slowly shrinks which maintains the high luminosity and causes a 

slow increase in the effective temperature of 3 x 10® K. The se- 
quence then cools rapidly and reaches equilibrium in about an hour. 

We extended our study to neutron stars with a radius of 20 km; 
bracketing the published work. In this case the reduced gravity 
should produce outbursts with lower peak luminosities and effective 
temperatures but the time scale for the outburst should be in- 
creased over that of the 10 kn evolution. Such was the case. 

For this evolution we chose an envelope mass of 2 x 10”^^ M. 

7 ® 

which gave a temperature and density at the CEI of 4 x 10 K and 

c .3 

4 x lO” gm cm , respectively. Because of the lower density, the 
evolution proceeds more slowly than for the 10 km case and it takes 

3 

nearly 10 sec for the peak temperature in the envelope to reach 

10 K. As before, this does not occur at the CEI, but a few zones 
closer to the surface. In addition, the lower gravity allows the 
envelope to expand to nearly 40 km during the evolution which pre- 
vents the temperature at the CEi from exceeding 4 x 10® K. At 

this temperature the ^^0(o,p)^^F and ^®0 (o,y)^'N 6 reactions are remov- 
ing catalytic nuclei from the CNO reaction sequence and are producing 
very little additional energy. This results in a most interesting 
phenomena. Once the shell source has reached inward to the CEI, 



the teitaeratures are high enough for a reasonable number of a- 

U IS 

reactions on 0 and 0 to occur. This is not true farther out 
in the envelope where the temperatures are lower. This results in 
the CNO reactions being able to cycle faster at lower temperatures 
causing an inversion in the hydrogen abundance. This can be seen 
in Figures 3 and 4 which show the tanperatures , and hydrogen abun- 
dances at two different times in the evolution. Therefore, some 
6 hours after the evolution began, the hydrogen abundance drops to 
zero about one- third of the way from the CEI to the surface. Be- 
cause helium is nearly depleted at the CEI the o-reactions on ^^0 
I*' 

and “0 also become less important and the principle source of 

energy comes from the 6^-decays followed by proton captures. Be- 
cause the CNO catalytic nuclei have their highest abundances at 
the edge of the region where hydrogen has become completely de- 
pleted, this is just the region where they are cycling the most 
rapidly and the burning front slowly moves inward and outward. How- 
ever, the hydrogen abundance is low and the runaway time scale 
for each zone is more rapid than the time scale for the inward 
diffusion of heat. This allows each zone to bum out before the 
next inner zone can flash, producing a very ragged curve for both 
the temperature versus time (Figure 5) and the luminosity versus 
time (Figure 6). If we had been able to use more zones, these 
curves would have been smoother. 

The peak rate of energy generation reached in this simulation 
15 -1 -1 

was only 10 erg gm s . The peak luminosity was 700 and the 

peak effective temperature barely exceeded 0.5 kev. It takes this 
sequence about or.» day to burn out all of the hydrogen in the shell 
source and for the luminosity to begin dropping. It then takes 
about 45 days to return to minimum. All of the hydrogen and helium 
are burnt to higher mass nuclei except for a thin shell of H of 
-13 

about 10 M which is still burning but on a much longer time 
SC3 Is ^ 

Because the ^^0(a,p)^^F and ^^0(a,Y)^^Ne rates are theoretical, 
we also investigated the effects of major changes in these rates by 
evolving one additional sequence with the identical initial condi- 
tions as in the previous sequence but with these rates set to zero. 

3 

It takes this simulation 2 x 10 sec for the peak temperature in 

8 

the envelope to reach 5 x 10 K. The rate of energy generation has 
.15 -1-1 

reached 10 erg gm sec . At this temperature a significant num- 

12 16 

ber of triple-cx and C(o,y) 0 reactions are occurring and the 
triple-i reaction is feeding new catalytic nuclei into the CNO. cycle. 
This increases the value of Z and counteracts the effects of the in- 
creasing number of nuclei being trapoed as ^^0 and ^^0. The tem- 

0 

cerature continues to increase finally reaching 10"K. The peak 
luminosity is 2.5 x 10^ and peak effective temperature is 1.2 kev. 


About 2 X '0^ sec later, hydrogen and helium burn out and the tem- 
perature starts to drop. It takes nearly 3 days for the effective 
te.mperature to fall below 0.1 kev. 

CONCLUSIONS 

We have found in this study that we can produce long time 
scale outbursts on neutron stars, if we assume low accretion rates 
and "cool" neutron stars. The time scales for these outbursts 

3 

range from 10 sec for the 10 km neutron star to about one day or 
longer for the 20 km neutron star. The peak temperatures and 
luminosities were inversely proportional to the radius of the neu- 
tron stars and our calculations (plus those noted earlier) suggest 
that the actual radii of most neutron stars must be closer to 10 
km than 20 km. On the other hand, the fast, soft. X-ray transients 
can be produced on larger radii neutron stars if such a wide range 
in neutron star radii is possible. 

We also produced flat topped outbursts similar to some of 
those observed. Such a theoretical outburst results when the ac- 
creted envelope has had enough time to reach thermal equilibrium 
before the outburst begins. 

Finally, we have been able to achieve outbursts in hydrogen 

rich material because the ^^0 (o,p)^”f and ^^0 (o,y)^^N 6 reactions 
act to remove catalytic nuclei from the CNO reaction cycle and, 
at high temperatures, restores the temperature dependence of the 
CNO reactions. In addition, at these same temperatures, the 
triple-a reaction is feeding new nuclei into the CN cycle which 
also keeps the rate of energy generation elevated over what one 
would predict if these reactions were not occurring. 
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